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We have considered a system of two superconducting charge qubits capacitively 
coupled to a microwave resonator. The dynamics of the Bloch vectors are inves- 
tigated for different regimes. By means of the Bloch vectors and cross dyadic we 
quantify the degree of entanglement contained in the generated entangled state. 
We consider different values of the system parameters to discuss the dynamics of 
the channel capacity between the qubits. We show that there is an important role 
played by initial state settings, coupling constant and the mean photon number on 
generating entangled state with high degree of entanglement and high capacity 
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1 Introduction 

Quantum entanglement is a quantum mechanical phenomenon in which the quantum 
states of two or more objects have to be described with reference to each other, 
even though the individual objects may be spatially separated. It is an essential 
component in many quantum information processing applications, such as quantum 
computation pQ, teleportation [3] and cryptography j2] and dense coding [U [5]. 
Therefore, it is essential to create and manipulate entangled states for quantum 
information tasks. The basic element of the quantum information is the quantum 
bit (qubit) which is considered as a two level system. Consequently, most of the 
research concentrate to generate entanglement between two level systems. Among 
these systems, superconducting charged qubits i.e. the Cooper pair box [6l [12] . 

In the last years the generated entangled states between Cooper pairs has been 
received a lot of attenuation. This is due to its properties as a two-level quantum 
system, which makes it a candidate as a qubit in a quantum computer [71 IB]. Also, 
recently these types of charged particles has been used to implement a Shor's factor- 
ization [9j, Deutsch-Jozsa algorithm [10] , the Grover search algorithm protocol [llj, 
and the quantum computing with a single molecular ensemble. On the other hand 
theses charged pairs have been used to realize a controlled phase gate [13] . All these 
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applications have encouraged different authors to study the entangled properties 
of these systems. As an example, the entanglement between two superconducting 
qubits has been generated by the interaction with nonclassical radiation [TJ]. Ro- 
drigues et al. have evaluated the entanglement of superconducting charge qubits by 
homodyne measurement [THIHE]. The entropy squeezing and the emission spectra 
of a single-cooper pair have been investigated [T7] . 

On the other hand, it is known that the Bloch vector gives one of the possible 
descriptions of N-level quantum states, because it is defined as a vector whose com- 
ponents are expectation values of some observables. So, we can use it to describe the 
density operator of any system [T5]. These vectors (called coherent vectors), play a 
central role in quantifying the degree of entanglement, where Englert [19] has used 
it to evaluate the entanglement dyadic which represent a measure of the degree of 
entanglement. Addressing this issue with Cooper pair box systems is the main aim 
of the present paper. 

In this article we propose to study a system of two-Cooper pairs interacting with 
a single cavity mode. We investigate the dynamics of the polarized vector for each 
individual qubit and the channel capacity of the entangled Cooper pair for different 
values of the system parameters. This paper is organized as follows. In section 2, 
we introduce our model and its solution. We devote section 3, to investigate the 
relation between the Bloch vectors and the quantum entanglement. In section 4, 
we investigate the effect of the structure of the initial state and the mean photon 
number on the transmission rate of information between two parties via the channel 
capacity. Finally, we end up with our conclusion. 

2 Model and its solution 

Our system consists of two superconducting charge qubits each coupled capacitively 
to a strapline resonator. Each charged qubit consists of a small superconducting is- 
land with a Cooper pair charge Q. This island connected by two identical Josephson 
junctions, with capacity Cj and Josephson coupling energy Ej, to a superconducting 
electrode. This system can be describe as a pair of two-level systems coupled to a 
simple harmonic oscillator. The charging energy of the qubits and their coupling to 
the resonator can be controlled by the application of magnetic and electric fields. 
By using the a rotating wave approximation, the system can be described by the 
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Hamiltonian |15j . 

1 2 
= hiu(a j a +7:) + {E 1 a z + E 2 t z ) + /l^ Ai(aa+ + aV - ), (1) 



i=l 



where a z = |e) 1 (e| — |5 , ) 1 ((y , | for the first qubit, = e)„(e| — |<7) 2 (#| for the second 
qubit, a) and a are the creation and annihilation operators of photons with frequency 
u, E\p are the charging energies of the qubits, and Ai : 2 are the resonator-qubit 
coupling terms. 

Our target of this section is to derive the time- dependent density matrix which 
enables us to discuss the statistical properties of the present model. For this reason 
let us assume that the initial state of the two charged qubits are prepared in a 
superposition state which can be written as r0 c (O)) = o|ee) + b \gg), \a\ 2 \ + \b\ 2 = 1, 
while the cavity field starts from a coherent state, |^(0)) = Y^Lol"]™) , 1 n = 
^f= exp(— ||a| 2 ). Now we can write the time evolution of the density operator in 
the form 

p(t)=W(t)p(0)Wt(t), (2) 

where lA(t) = exp (—iHt/h) is the time-dependent unitary operator and p(0) = 
|^(0)) (^(O) I . Since the invariant sub-space of the global system can be considered as 
a set of complete basis of the atom-field. Since, we are interested in the contribution 
of the Bloch vector in quantifying the degree of entanglement, we express the density 
operator p c (t) of the system charged system by means of the Bloch vector for each 
qubit and the cross dyadic. After tracing out the state of the field one gets the state 
of the charged qubit as, 

p c {t) = -(1 + s ■ a 1 + t ■ t 1 + a -%rrowC ■ t 1 ), (3) 

where a = (a x ,a y ,a z ), = (r z , r y , r z )are the Pauli matrices for the first and 
second qubit respectively, while s = (s x ,s y ,s z ), t = (t x ,t y ,t z ) and are the Bloch 
vectors for the first and the second qubits and l @rrowC is the cross dyadic, are 
defined as follows: 

s x (t) = A n+1 (ty n (t) + B n+1 (t)D* n (t) + C n (t)A* n+1 (t) + D n (t)B* n+1 (t), 

s z (t) = \A n (t)\ 2 + \B n (t)\ 2 -\C n (t)\ 2 -\D n (t)\ 2 , (4) 
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while for the second qubit, 



t x (t) = A n+1 (t)B* n (t) + B n (t)A* n+1 (t) + C n+1 (t)D* n (t) + D n (t)C* n+1 (t), 
t y (t) = i [-A n+1 (t)B* n (t) + B n (t)A* n+1 (t) - C n+l (t)Dm + D n (t)C* n+1 (t)} , 
t,(t) = \A n (t)\ 2 -\B n (t)\ 2 + \C n (t)\ 2 -\D n (t)\ 2 , (5) 

In addition to the elements of the cross dyadic, l @rrowC 



= i 
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C-zz 



A n+2 (t)D* n (t) + B n (t)C* n (t) + C n (t)B* n (t) + D n (t)A* n+2 (t), 
i [-A n+2 (t)D* n (t) + B n (t)C* n (t) - C n (t)B* n (t) + D n (t)A* n+2 (t)} , 
A n+1 (t)C* n (t) - B n+1 (t)D* n (t) + C n (t)A* n+1 (t) - D n (t)B* n+1 (t) 
i [-A n+2 (t)D* n (t) - B n (t)C* n (t) + C n (t)B* n (t) + D n (t)A* n+2 (t)] , 
-A n+2 (t)D* n (t) + B n (t)C*(t) + C n (t)B*(t) - D n (t)A* n+2 (t), 
i [-A n+1 (t)C* n (t) + B n+1 (t)D* n (t) + C n (t)A* n+l (t) - D n (t)B* n+1 (t)} 
A n+1 (t)B* n (t) + B n (t)A* n+1 (t) - C n+1 (t)D* n (t) - D n (t)C* n+1 (t), 
i [-A n+l {t)B* n {t) + B n {t)A* n+1 {t) - C n+1 (t)D* n (t) + D n (t)C* n+1 (t)} 

\A n (t)\ 2 -\Bn{t)\ 2 -\C n {t)\ 2 + \Dn{t)\\ 
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and, 7n = y/l + n, f3 n = y/n + 2 fi n = \(5 n + y/ 5 2 n -4 A 2 ), i/ n = \(5 n - 
- 4A 2 ) tf B = (1 + i? 2 )( 7 2 + /3 2 ) , A n = (1 - i? 2 ) 2 /^ 7 2 and i2 = £ 

3 Dynamics of Bloch Vectors 

In this section, taking into account the possible extension of the methods which 
will be discussed in this work to higher-dimensional situations and to mixed states, 
the present strategy will certainly need to be modified and integrated by developing 
more sophisticated tools, including the use of multiple operations to yield properly 
defined entanglement and generalizations to mixed-states cases. At this end, we can 
say that, the qubit |^} = a|0) + can be represented as a point (9, 0) on a unit 
sphere called the Bloch sphere. Define the angles 9 and by letting a = cos(|) and 
b = e 1 ^ sin(|), then |-0) is represented by the unit vector (cos sin 9, sin sin 9, cos 9) 
called Bloch vector. 

Here, we investigate the dynamics of the Bloch vectors (Polarized vector), for 
different values of the mean photon number n, the relative coupling R and different 
classes of the initial state of the charged system. In Fig. (1), we assume that the 
charged qubits are prepared in the excited state p e (0) = |ee)(ee| and the mean 
photon number n = 20, while the relative coupling R = A2/A1 takes different values. 
In Fig. (la), we consider a small value of R, which means that the coupling between 
the second qubit and the cavity field is very weak. It is clear that, the Bloch vector 
for the second qubit is bigger than the Bloch vector for the first qubit. This is due 
to the weak coupling between the second qubit and the cavity field. So, the chance 
that the second qubit interacts with the field is very small. On the other hand, as 
soon as the interaction goes on, we can see that the Bloch vector for the first qubit 
decreases rapidly. This is due to the strong coupling between the first charged qubit 
and the cavity mode. So for this situation, there is an entangled state generated 
between the first qubit and the filed, while the second qubit is almost factorized. 
This phenomenon is not observed in Fig. (lb), where we consider large values of 
the coupling R. For the numerical values we consider R = 0.9, i.e Ai ~ A2. In 
this case the two Bloch vectors increase and decrease together. This means that the 
probability that both qubits interact with the field is almost the same. Also, we can 
notice that, this choice causes a shrink of the Bloch vectors, where the maximum 
values almost 0.6. 

It is interesting to mentioned to the fact that, the maximum entangled state as 
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Figure 1: The Bloch vectors as functions of the scaled time. The two-qubit system 
is prepared initially in an excited state \ei, €2), and the field starts from a coherent 
state with a mean-photon number n = 20. The solid curves for the first qubit while 
the dot curves for the second qubit. The relative coupling is (a) for R = 0.003 (b) 



Bell states has zero Bloch vectors. So, as the Bloch vectors decrease, the possibility 
for obtaining an entangled state with high degree of entanglement increases. So, by 
using a suitable values of the coupling constant, one can control the Bloch vectors 
and consequently it is possible to generate an entangled state with high degree of 
entanglement. In Fig. (2), we plot the Bloch sphere, with radius equal 0.4 units, at 
specific time in the interval [10.2 — 11.6], where we consider that the charged qubits 
are prepared in the excited state and R = 0.9. The amplitude and the direction of 
the Bloch vector for the first qubit are shown. One can see that at t = 10.2, the 
Bloch vector parallel to y-axis and it has a large value. As the time increases (say 
t = 10.3, 10.4), the Bloch vector inside the sphere and it could great or equal 0.4. 
As the time increases further, the Bloch vector appears in a different direction and 
its amplitude shrinks, this is shown for t = 10.5. As time increases more the length 
of the Bloch vector increases and its direction change. 

In Fig. (3), we assume that the charged system is prepared in the partially 
entangled state a|ee) + b\gg^. For small values for the relative coupling the general 
behavior of the two vectors is somehow similar with small differences. Comparing 
Figs. (la) and (2a), we see that the minimum point of the Bloch vector is decreased 
for the partial entangled state, which means that, using this class of the initial state, 
it will be more efficient to generate entangled state with high degree of entanglement. 
For large values of the relative coupling constant, R = 0.9, the behavior of the two 



R = 0.9. 
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Figure 2: The Bloch sphere for a charged system prepared initially in an excited 
state with n = 20, R = 0.9. The thick line represents the length of the Bloch vector 
for the first qubits. 

Bloch vectors are similar. For a specific time one of the Bloch vectors reaches to 
zero, in this case there is a maximum entangled state generated between the two 
charged qubit, 

In Fig. (3c), we consider the mean photon number n = 10 and R is small. In 
this case, we can see that the minimum point of the Bloch vectors is very small 
compared with that depicted in Fig. (3a). Also, the Bloch vector for the second 
qubit is shrunk more for small values of the mean photon number. This phenomena 
is clearly appeared in Fig. (3d), where we consider a large value of the coupling 
constant. In additional to the coincidence of the behavior of the Bloch vector for 
the two qubits, the two vectors are shrinking more. So, in this case the possibility 
of generating entangled state with high degree of entanglement is increased. 



4 Degree of Entanglement 

To quantify the degree of entanglement between the two charged qubits, we use a 
measure which is defined by means of the Bloch vectors and the cross dyadic. We 
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define the entangled dyadic as 



) rrowE 



\rrowC — s t 



where l @rrowc is a 3 x 3 matrix which is defined by(jHJ) and t is also 3x3 matrix 
whose elements can be obtained from (j3j) and (jSJ). The degree of entanglement is 
defined by 

DoE = tr{ l @rrowE T - L @rrowE} (9) 

where ©rrowE is the transpose matrix of the dyadic ©rrowE 

The amount of entanglement between the entangled charged qubits is shown in 
Fig. 4, in which we consider different regimes. In the first regime, we assume that 
the two qubits are prepared initially in the following excited state IV'(O)) = |ee) 
while the field starts from a coherent state with a mean photon number n = 20. In 
Fig. (4a), we see that for a small value of the coupling constant (say R = 0.003), 
there is no any quantum correlation between the two charged qubit expect on the 
interval [19.2-40.5], where in this interval the Bloch of the second qubit decreases. 
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Figure 3: The same as Fig.(l) but the charged qubits are prepared initially a partial 
entangled state. 
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Figure 4: The degree of entanglement for a charged system prepared initially in an 
excited state with n = 20. (a) For R = 0.003 (a) excited (b) R = 0.9. 



This means that in this intervals, the three subsystems (the two charged qubits 
and the field) interact with each others. Also, from this figure, we can see that the 
maximum amount of entanglement is obtained at the minimum point of the Bloch 
vectors of both qubits. This phenomena, also, is shown in Fig. (46), where we 
consider R = 0.9. As soon as the interaction starts, (scaled time is greater than 
zero), the entangled state starts to be generated. It is obvious to realize that the 
development of the entanglement depends on the dynamics of the Bloch vectors. 

The effect of different setting of the initial state of the charged qubit is seen in 
Fig. (5), where we assume that the system is prepared initially in a partially entangled 
state. Fig. (5a), is devoted to consider a weak correlation case (say R = 0.003). 




10 20 30 40 50 60 10 20 30 40 50 60 
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Figure 5: The degree of entanglement for a charged system prepared initially in 
partially entangled states with n = 20. (a) For R = 0.003 (a) excited (b) R = 0.9. 
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From this figure it is clear that the degree of entanglement decreases with time until 
it reaches a minimum point at Ait ~ 15. From Fig. (3a), we see that the Bloch 
vector for both qubit are maximum and consequently the degree of entanglement 
will be minimum according to the definition (Q. In the interval [40 — 60], the degree 
of entanglement starts to increase another time, where in this interval of time the 
Bloch vector decrease (see Fig. (3a)). For a strong correlation the dynamics is seen 
in Fig. (56). It is clear that the degree of entanglement is much stronger, where 
the efficiency of generating entangled state in this case is much higher. Also, the 
evolvement of the degree of entanglement depends on the evolvement of the Bloch 
vector for each qubit. 

In this context it is very important to investigate the evolvement of the degree 
of entanglement for different values of the mean photon number n. We choose the 
case for the weak coupling, where R = 0.003. In Fig. (6a), we assume that the 
system has been processed in advance in an excites state. The usual effect of the 
mean photon number is seen where the Rabi oscillation is shifted to the left. In 
this case the maximum amount of entanglement is smaller than that depicted for 
n = 20 (Fig. (4a)). But on the other hand, there is an entangled state is generated 
in the interval [50 — 60], while on the corresponding interval the charged system 
behaves as a product state. In Fig. (66), we plot the degree of entanglement for 
a charged qubit is prepared in partially entangled state. In this case the degree of 
entanglement is much better than that shown in Fig. (5a), where n = 20. So, For 
the weak interaction one can generate an entangled state between the two charged 
qubits by reducing the number of photons in the cavity mode. 

5 The channel Capacity 

In this section we investigate the effect of the structure of the initial state and the 
mean photon number on the transmission rate of information from a sender (Alice), 
to a receiver (Bob). This task can be performed by employing the dense coding 
protocol [3], [2QI The main idea of this protocol is that, Alice and Bob share 
an entangled qubit pair. They used it as a channel, where Alice can encode two 
classical bits in her qubit by using local operation. Alice sends her qubit to Bob, 
who will try to decode the information. The amount of information gained by Bob 
depends on the capacity of the channel. In this context, we try to show how the 
capacity of the channel and hence the rate of transmit information depend on the 
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Figure 6: The degree of entanglement for a charged system prepared initially in 
partially entangled states with n = 10, R = 0.003. For (a) excited (b) partially 
entangled state. 

structure of the initial states, the relative coupling and the mean photon number. 
For N x M state system, the dense coding capacity is given by 

C = logD A + S(p B )-S(p AB ), (10) 

where Da = N, p B = ^^{pab} an d <S{-) is the Von Numann entropy. Since we 
consider entangled two qubits, then our system is in 2 x 2 dimension. Fig. 5, shows 
the behavior of the capacity of the entangled quantum channel, where we consider 
the two charged qubits are in excited state. We investigate the effect of the coupling 
constant, the mean photon number and the setting of the initial state of the charged 
system. Fig. (7a), shows the behavior of the channel capacity for a charged system 
with a strong coupling with the field, where we assume that R = 0.9 and the mean 
photon number n — 10. From this figure it is clear that the channel capacity for a 
charged system is prepared initially in a partially entangled state is much better than 
that if we start with a charged qubit in an excited state (product state). The effect 
of the mean photon number is seen if we take a look at Fig. (7b), where we consider 
n = 10. From this figure, we can see that the Rabi oscillations shrunk and shifted 
to the right as time increases. Due to the shrunk of the oscillations, we can see that 
the channel capacity increases a little bit as n increases. From Fig. (46) and Fig. (76), 
we can see that there is a strong relation between the degree of entanglement and 
the channel capacity and hence the rate of transform data. This maximum value of 
the channel capacity is obtained at a large degree of entanglement, as an example in 
the interval [0.2,0.8]. Also, from equation (fTUI) . the capacity of the channel does not 



11 



only depend on the channel p atn but also on the individual state for the single qubit 
p b . It is clear on the interval [10-20] as an example, the degree of entanglement is 
almost zero, but the capacity does not vanish. 



6 Conclusion 

In this paper we have studied in a non-standard way the dynamics of the Bloch 
vectors for two charged qubits. We investigate the dynamical behavior of the Bloch 
vector for each qubit. The shrunk, extension and the direction of these vectors are 
examined for different parameters of the charged system and the cavity field. We 
show for strong coupling between the filed and the charged system, the Bloch vector 
for the two qubits have the same behavior. But if one qubit has a weak coupling 
with the field (as the second qubit in our case), the behavior of the two Bloch vectors 
is different. 

Using the density matrix technique, we predict the existence of entangled states, 
where we consider an entanglement measure which depends on the two Bloch vectors 
and the cross dyadic and is called entangled dyadic. The relation between the 
evolvement of the Bloch vector and the degree of entanglement is shown, where for 
large values of the Bloch vectors, the degree of entanglement is minimum. It is 
shown that, for a charged qubit prepared initially in an entangled state, the amount 
of entanglement is much larger than that for any other choices. The role played 
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Figure 7: The capacity of the channel where R — 0.9 The solid and dot curves 
for a charged system is prepared initially in a partially entangled and excited states 
respectively, (a) for n = 10, and (b) For n = 20. 
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by the coupling constant and the mean photon number is cleared for generating 
entangled states and improving its degree of entanglement. Decreasing the mean 
photon number of the cavity mode is important to generate entangled state even 
the coupling between the qubits and field is weak. 

The behavior of the channel capacity between the two charged qubits is ex- 
amined, using different regimes of preparing the initial state of the system. It is 
shown that the possibility of generating entangled state with high capacity and 
consequently high rate of transmission of information is much better if we start 
with a partially entangled state for the system. Also, the mean photon number 
plays a central role in the efficient of the channel capacity, where for small values of 
the mean photon numbers, the channel capacity is increased and consequently the 
transmission rate of information. 
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